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> ' ABSTRACT 

! We present a spectroscopic analysis of five stellar streams ('A', 'B', 'Cr', 'Cp' and 'D') 

. as well as the extended star cluster, EC4, which lies within Stream 'C, all discovered 

in the halo of M31 from our CFHT/MegaCam survey. These spectroscopic results 
were initially serendipitous, making use of our existing observations from the DEcp 
Imaging Multi-Object Spectrograph mounted on the Keck II telescope, and thereby 
' emphasizing the ubiquity of tidal streams that account for ~70% of the M31 halo stars 

in the targeted fields. Subsequent spectroscopy was then procured in Stream 'C and 
Stream 'D' to trace the velocity gradient along the streams. Nine metal-rich ([Fe/H]~- 
0.7) stars at t>h c i = —349.5 km/s, cr ViCorr ~ 5.1 ± 2.5km/s are proposed as a serendipi- 
tous detection of Stream 'Cr', with followup kinematic identification at a further point 
along the stream. Six metal-poor ([Fe/H]~-1.3) stars confined to a narrow, 15 km/s 
velocity bin centered at t>hd — —285.6 km/s, a ViCorr = 4. 3^4 km/s represent a kine- 
matic detection of Stream 'Cp', again with followup kinematic identification further 
along the stream. For the cluster EC4, candidate member stars with average [Fe/H]~- 
1.4 ([Fe/H] spec =-1.6), are found at Vhd — —285 km/s suggesting it could be related 
to Stream 'Cp'. No similarly obvious cold kinematic candidate is found for Stream 'D', 
although candidates are proposed in both of two spectroscopic pointings along the 
stream (both at ~ — 400km/s). Spectroscopy near the edge of Stream 'B' suggests a 
likely kinematic detection at t>hoi ~ —330 km/s, (J VtCorr ~ 6.9 km/s, while a candi- 
date kinematic detection of Stream 'A' is found (plausibly associated to M33 rather 
than M31) with v^el ~ —170 km/s, o v , corr — 12.5km/s. The low dispersion of the 
streams in kinematics, physical thickness, and metallicity makes it hard to reconcile 
with a scenario whereby these stream structures as an ensemble are related to the 
giant southern stream. We conclude that the M31 stellar halo is largely made up of 
multiple kinematically cold streams. 

Key words: 



1 INTRODUCTION 



Observatory, which is operated as a scientific partnership among 
The data presented herein were obtained at the W.M. Keck the California Institute of Technology, the University of Califor- 
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Stellar streams represent the visible debris of small galax- 
ies being cannibalized by large galaxies, memorials to the 
merging process by which the halos of galaxies are built 
up. The best known examples of streams in the Milky Way 
(MW) have recently been mapped far more extensively by 
the SDSS-DR5 by Belokurov et al. (2006, 2007): the tidally 
stripped stars and globular clusters associated with the 
Sagittarius dwarf spheroidal and the Low Latitude stream, 
along with a newly discovered "Orphan Stream" so named 
for its lack of obvious progenitor. In M31, thanks to our abil- 
ity to efficiently map vast regions of the halo, the number 
of discovered giant streams already outnumbers that of the 
MW (Ibata et al. 2007). However, of the eight stellar streams 
that have been identified in the halo of M31, only one has 
plausibly been identified to a dwarf satellite: the loop con- 
necting to NGC 205 (McConnachie et al. 2005). This sug- 
gests that the other streams could represent an additional 
seven 'uncataloged' satellites, although some of the streams 
might be produced by a common progenitor, as suggested 
by models of Fardal et al. (2007, 2008) for the M31 Giant 
Southern Stream, or in a similar fashion to the Sgr dwarf 
and its numerous wraps around the Milky Way. It is impor- 
tant to characterize their orbits, metallicities and masses, to 
understand what their progenitors must have been. 

The existence of stellar streams tells us that a progen- 
itor galaxy has undergone significant mass loss. This is due 
to a combination of its orbit and its phase of evolution - the 
amount of dark matter mass the galaxy has lost (Penarrubia 
et al. 2008a, b) . Satellites on circular orbits are harder to dis- 
rupt, but if they are massive enough (i.e. of the order of the 
LMC), dynamical friction will bring them close to the host 
galaxy centre, where the interactions with the disc will lead 
to their tidal disruption - e.g., Penarrubia et al. 2007). In 
addition, to form a stream one has to remove most of the 
dark matter halo (~90-99%). The most important param- 
eter that controls the mass loss rate of a dSph is the peri- 
centre distance (the orbital eccentricity is of second order). 
On the other hand dwarfs with highly elliptical orbits spend 
a lot of time near apocentre where they are unlikely to be 
disrupted by the host. It is therefore not immediately obvi- 
ous that streams represent preferred types of orbits on av- 
erage. However, a number of theoretical studies have shown 
that significant information about the orbital properties of 
the progenitor galaxy can be derived from the streams (e.g, 
Ghigna et al. 1998; Helmi et al. 1999a). 

Streams can be much more informative to study than 
dwarfs because their orbits can be directly traced and con- 
strained. Fellhauer et al. (2006) were able to accurately con- 
strain the shape of the Galactic potential through the bi- 
furcation of Sagittarius streams in Belokurov et al. (2006). 
Understanding the range of orbits of satellites to large galax- 
ies will help us to understand how the halos of these galaxies 
formed. This is especially interesting in light of M31's huge 
stellar halo reflecting the dark matter dominated halo out 
to > 150 kpc (e.g., Irwin et al. 2005, Gilbert et al. 2006, 
Ibata et al. 2007). However streams can also be much harder 
to analyse observationally: the distances are problematic, 
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there's a much lower spatial density and they have a larger 
extent so that observational sampling is not trivial. There is 
also the difficulty to infer the membership of different stream 
pieces, especially if we expect different chemical signatures 
due to metallicity gradients in the progenitor system (e.g., 
Ibata et al. 2007). 

Streamy/blobby structures are individually interesting 
and constraining for the halo formation. They can represent 
the only traceable product of long disintegrated progenitors, 
yet retain a coherent body for statistical analysis. Streams 
can provide important clues on the structure of the pro- 
genitors (e.g. metallicity gradients, mass to light - M/L) 
as well as on the shape of the host dark matter halo (e.g. 
prolate versus oblate) (Martinez-Delgado et al. 2008). Fu- 
ture study of these structures will be able to put them in a 
much better "near-field cosmology" context, eventually un- 
derstanding their ages and chemical histories. However it is 
important to uncover and study them now, even in limited 
capacities necessitated by the small numbers of spectroscop- 
ically identifiable stars and HR-diagram depths, so we can 
build our models on the most complete context. 

We have initiated a spectroscopic survey of the new 
streams found in M31's halo using the DEep Imaging Multi- 
Object Spectrograph on Keck II to derive radial velocities 
and metallicities of red giant branch (RGB) stars. In this 
contribution, we discuss spectroscopic pointings in each of 
streams 'C' and 'D' which we obtained by serendipity, since 
the spectroscopy was taken prior to knowledge of the pho- 
tometrically discovered streams, as well as followup spectro- 
scopic pointings in both of these streams. We also analyze 
spectroscopic data from Koch et al. (2008) lying within the 
Ibata et al. (2007) streams 'A' and 'B'. 



2 OBSERVATIONS AND ANALYSIS 

The spectroscopic fields along the M31 minor axis discussed 
in this paper are highlighted in Figure 1, two lying at 
~35 kpc, three lying at ~60 kpc, one at ~80 kpc, and one 
at ~120 kpc projected from the center of M31. The fields 
cover the four stellar streams presented in the Ibata et al. 
(2007) M31 extended halo analysis, called Stream 'D', 'C\ 
'B', and 'A' respectively. In these M31 halo images, it can 
be seen that Stream 'C' has significantly different morphol- 
ogy as a function of metallicity, a more metal-rich compo- 
nent dominating the structure, with a more irregular shaped 
metal-poor component to the east. While it was not pro- 
posed initially, our evidence in this work suggests the two 
structures may be distinct systems, and we refer to these as 
two separate streams, Stream 'Cr' for the metal-rich com- 
ponent and Stream 'Cp' for the metal-poor component (this 
issue is explored in detail in § 3.3 and Figs. 10 & 11). 

Multi-object spectroscopic observations with the Keck- 
II telescope and the DEep Imaging Multi-Object Spectro- 
graph - DEIMOS (Davis et al. 2003) were obtained in pho- 
tometric conditions with ~ 0.8" seeing in Sept. 2004 and 
2005. Target stars were chosen by colour/magnitude selec- 
tion as described in Ibata et al. (2005), first selecting likely 
RGB stars in M31 over all metallicities, and filling space 
with any other stellar objects in the field. Two spectro- 
scopic masks (F25 and F26 from the table in Chapman 
et al. 2006) targeted the field of an extended cluster EC4 
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Figure 1. The locations of the spectroscopic measurements (red circles) are overlaid on the imaging data from the INT and CFHT 
telescopes presented in Ibata et al. (2007). Streams 'B', 'C and 'D' intersect the minor axis of M31 approximately perpendicularly at 
~ 80 kpc, ~ 60 kpc and ~ 35 kpc, respectively. The irregular turquoise line demarks the CFHT survey region. The scale of the diagram 
is shown by the circle segment (dashed line - marking a projected radius of 100 kpc), as well as by the ellipse segment (continuous line - 
showing a 50 kpc ellipse of axis ratio 0.6). Two of the newly discovered dwarf spheroidals from our survey are visibile in the same region, 
And XII (Martin et al. 2006, Chapman et al. 2007), and And XV (Ibata et al. 2007, Letarte et al. 2008). 



(Mackey et al. 2006), which were found after the fact to 
be spanning Stream 'Cr' and Stream 'Cp'. A combined to- 
tal of 212 independent stars in both masks were observed 
in standard DEIMOS slit-mask mode (Davis et al. 2003) 
using the high resolution 1200 line/mm grating, and 1" 
width slitlets. Ten of these target stars were specifically se- 
lected from HST photometry of EC4 to lie within the clus- 
ter. Our instrumental setting covered the observed wave- 
length range from ~ 0.70-0.98 /im. Exposure time was 60 
min, split into 20-min integrations. The DEIMOS-DEEP2 
pipeline (Newman et al. 2004) designed to reduce data of this 
type accomplishes tasks of debiassing, flat-fielding, extract- 
ing, wavelength-calibrating and sky-subtracting the spectra. 
The same settings were used to target a halo field which was 
found after the fact to lie in Stream 'D'. 89 stars were ob- 
served in this mask (F7). 

The Ibata et al. (2007) imaging discovery of the new 
streams, coupled with obvious kinematic detection of the 
Stream 'Cr' and Stream 'Cp' in our existing spectroscopic 
observations (described in subsequent sections) prompted 
the followup study of these structures. On October 8, 2007, 
additional DEIMOS masks were obtained further along the 
Stream 'C and Stream 'D', as identified in Fig. 1 and Ta- 
ble 1. These observations were obtained under ~ 1 arcsec 
seeing, and cloudy conditions. For the field F36, we obtained 



3x20min integrations on a mask with 113 targeted stars, 
while field F37 was observed for 4 x 20 min integrations with 
138 targeted stars. These observations used the lower reso- 
lution 6001/mm grating to achieve a higher S/N in the con- 
tinuum for the fainter stars, and resulting in a resolution of 

o 

~3 A estimated from the width of sky lines. 

The radial velocities of the stars in all these fields were 
then measured with respect to spectra of standard stars ob- 
served during the observing runs. By fitting the peak of the 
cross-correlation function, an estimate of the radial velocity 
accuracy was obtained for each radial velocity measurement. 
The accuracy of these data, as estimated from the Calcium 
Triplet (CaT) cross-correlation, varies with magnitude, hav- 
ing uncertainties of < lOkrns -1 for most of the stars. The 
CMDs, velocity errors, velocity histograms, and metallicities 
for these fields are shown in Figures 2 & 3 for Stream 'C, 
and Figures 4 & 5 for Stream 'D'. Spectroscopic metallici- 
ties quoted in these tables are calculated from the equivalent 
widths of the Call triplet lines, as described in Ibata et al. 
(2005). 

Finally, the two furthest streams along the minor axis 
from Ibata et al. (2007), Stream 'A' and Stream 'B' ly- 
ing at 120 kpc and 80 kpc respectively, have serendipitous 
DEIMOS spectroscopic pointings lying in their edge regions 
from Gilbert et al. (2006) and Koch et al. (2008) (fields M8 
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and M6 respectively). Reduction and analysis of these two 
fields is detailed in Koch et al. (2008) . Figures 6 & 7 show the 
CMDs and velocity/metallicity distributions for Stream 'B' 
& Stream 'A' respectively. 

We address Galactic contamination to our spectroscop- 
ically identified stars in a manner identical to Koch et al. 
(2008), using a combination of V — I radial velocity and 
the equivalent width (EW) of NalA8i83,8i95 which is sen- 
sitive to surface gravity, and is accordingly very weak in 
M31 RGB star spectra, but can be strong in Galactic dwarfs 
(Schiavon et al. 1997). At velocities v hc i < —150 km/s, very 
few RGB candidates show any significant NalA8i83,8i95 ab- 
sorption lines, whereas stars with —150 to km/s veloc- 
ity show strong Nal absorption on average, consistent with 
the findings of Guhathakurta et al. (2006), Chapman et al. 
(2006), Gilbert et al. (2006) and Koch et al. (2008). For 
this study, we impose the additional constraint of remov- 
ing all stars from the halo sample with Vhoi > —150 km/s, 
and we remove any stars from our sample which have 
a summed EW(NalA8i83,8i95)> 0.8 in the velocity range 
Vhci < —150 km/s. 

The properties of all candidate M31 halo (and stream) 
stars in these fields are listed in Tables 2-7, including coor- 
dinates, velocities, spectroscopic and photometric metallici- 
ties, and V,I photometry. 

2.1 Observations of the cluster, EC4 

The two minor-axis stellar halo fields, F25/F26, lying 
serendipitously in Stream 'C' were observed with the addi- 
tional goal of constraining the kinematics of an "extended 
star cluster", EC4. These extended, luminous objects in the 
outskirts of M31 represent a population with ~ 15 - 60 kpc 
projected radii, large half-light radii for GCs with luminosi- 
ties near the peak of the GC luminosity function (Huxor et 
al. 2005, 2008). As such, they are dissimilar to any other 
known clusters in the Milky Way or M31, and begin to fill 
in the gap in parameter space between classical globular 
clusters and dwarf spheroidal galaxies. The "faint fuzzies" 
discovered in NGC 1023 (Larsen & Brodie 2000, Brodie & 
Larsen 2002), and the similarly diffuse objects in the ACS 
Virgo Cluster survey (Peng et al. 2006), may represent a sim- 
ilar class of cluster. EC4 was discovered within the CFHT- 
MegaCam survey (Huxor et al. 2008) at a projected radius 
of 60 kpc (coordinates in Table 1). 

Ten stars were selected to lie within EC4 by their CMD 
colours. Of these ten, one was subsequently identified as a 
galaxy in the HST/ACS image (Mackey et al. 2006), another 
lies well off the cluster RGB, likely due to contamination 
in the ground-based CFHTmegacam imaging affecting the 
colour measurements, while a third lies at ~5 core radii and 
has a low probability of being associated to the cluster. The 
remaining seven stars are candidate EC4 members, lying 
within 3 core radii of the cluster centre and falling along the 
top of the narrow RGB from the HST imagery (Collins et 
al. in preparation). 

In Fig. 2, the MegaCam photometry is shown for EC4 
stars for consistency with the stream data. However com- 
parison with the very narrow RGB from HST photometry 
in Mackey et al. (2006) shows that crowding affects the 
ground-based accuracy, since the scatter on the Megacam 
CMD is much larger than the difference in photometric er- 



rors (HST to Megacam) would warrant. In analyzing the 
DEIMOS spectroscopy, one targeted star lay exactly on the 
EC4 CMD, but had a discrepant velocity from the others. 
Closer examination of the spectrum revealed good detec- 
tions of the first and second CaT lines with a velocity of 
-277.3km/s, whereas the automated software pipeline de- 
rived a cross-correlation fit to larger skyline residuals. We 
include this star in the catalog as a viable member of EC4. 

2.2 Velocity accuracy: repeat measurements of 
stars in fields F25 & F26 

While CaT fitting errors suggest relatively small velocity 
errors, an independent check can be made on the 56 ra- 
dial velocities of stars lying in spectroscopic masks of both 
fields F25 and F26. Velocity differences are shown in Fig. 8 
highlighting those corresponding to Stream 'O (3 stars), 
Stream 'D' (2 stars), EC4 (3 stars), background halo (2 
stars), and Milky Way foreground (46 stars). A systematic 
shift from night 1 to night 2 of 3 km/s was found over all 
velocities, and this has been removed as a constant. Agree- 
ment between observing nights for these stars is then gen- 
erally found within the la errors of the radial velocity mea- 
surements, suggesting that no significant skew from mask 
misalignments or systematic errors are present from instru- 
mental setup night to night. The dispersion in velocity dif- 
ferences is ~6 km/s for both the M31 sample and the Milky 
Way sample taken separately, which is comparable to the 
typical velocity measurement error of an individual star. For 
these 56 stars, we have taken as the radial velocity the error- 
weighted average of the measurement from the two nights. 



3 RESULTS 

3.1 Global kinematics and metallicities of the 
fields 

Figures 2-7 (left panels) show the CFHT-megacam colour- 
magnitude diagrams (CMD) and radial velocity uncertain- 
ties of the observed stars in the Stream 'C, 'D', 'B', and 'A' 
fields, while the right panels show velocity histograms and 
photometrically derived metallicities, [Fe/H] = \og(Z/Zo), 
are computed for the stars by interpolating between 10 Gyr 
old Padova isochrones (Girardi et al. 2004). The fixed age 
adopted for metallicity comparison to the fiducials will of 
course introduce a systematic uncertainty if they aren't all 
old. Given that younger populations have been detected in 
the halo of M31, and that these streams could represent 
progenitors with a range of properties, we provide an es- 
timate of the age variations on the [Fe/H] determinations. 
If 5 Gyr old isochrones are used, the apparent metallicity 
would shift by 0.2 dex more metal-rich. The average distance 
modulus of M31 is adopted, 24.47, 785 kpc (McConnachie 
et al. 2005). Stars which are unlikely to be contaminated by 
foreground Milky Way (vhel < —150 km/s) are highlighted. 
While spectroscopic metallicities are quoted in the Tables, 
we do not use them for analysis here as the errors on in- 
dividual star measurements are so large as to broaden the 
typical [Fe/H] distribution by a factor of three for a stream 
kinematic structure. While the photometric [Fe/H] determi- 
nations are highly model dependent, the distributions for a 
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Figure 2. Left panel: CFHT-MegaCam color-magnitude diagram and radial velocity uncertainties of the observed stars in the 
Stream 'Cr'/Stream 'Cp' fields. Stars likely belonging to Stream 'Cr' (red), Stream 'Cp' (cyan), and EC4 (blue) are highlighted. The 
bright and very blue star belonging to Stream 'Cp' is unusual for M31 and its properties are described in the text. The fiducial RGBs 
correspond to, from left to right, NGC 6397, NGC 1851, 47 Tuc, NGC 6553 which have metallicity of [Fe/H]=-1.91, -1.29, -0.71, and 
-0.2, respectively. These fiducials have been shifted to the average distance modulus of EC4, 24.47 (785 kpc - Mackey et al. 2006). Right 
panel: The velocities of observed stars in the F25/F26 fields are shown as a histogram, with EC4 member stars highlighted as a heavy 
histogram. The stellar halo velocity dispersion (<r„=125 km/s) from Chapman et al. (2006) is shown normalized to the expected 9 halo 
stars at this position from Ibata et al. (2007). To differentiate EC4 stars from the field, we additionally plot the velocities against their 
radius from the EC4 center (Table 1), referencing the symbols to the CMD plot. Photometrically derived [Fc/H] is shown as a function 
of radial velocity, again referenced in symbol type to the CMD plot. 
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Figure 3. The same as for Fig. 2 for a field further along the Stream 'C structure (F36). No similarly obvious kinematic peaks are 
detected as in Fig. 2. In this Figure only we highlight stars potentially belonging to Stream 'Cr' (red) and Stream 'Cp'(cyan) entirely 
by their velocity ranges consistent with the spikes found in Fig. 2. Clusters of metal-poor and metal-rich stars are indeed found in 
these velocity ranges consistent with belonging to stream'Cp' and stream'Cr' respectively, although it is difficult to separate them from 
field/halo stars. 
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Figure 4. left panel: CFHT-megacam color-magnitude diagram and radial velocity uncertainties of the observed stars in the Stream 'D' 
field (F7). The fiducial RGBs are as in Fig. 2, shifted to the same average distance modulus of M31, 24.47, 785 kpc (McConnachie et 
al. 2005). Large symbols represent stars unlikely to be contaminated by foreground Milky Way (v^ e i < —150 km/s). Symbols are 
further highlighted which have photometric mctallicities consistent with the [Fe/H] distribution of Stream 'D' in Ibata et al. (2007): 
— 1.7 < [Fe/H] < —0.7. right panel: The velocities of observed stars in the Stream 'D' field are shown as a histogram, with possible 
Stream 'D' stars from the left panel highlighted as a filled histogram. Photometrically derived [Fe/H] is shown as a function of radial 
velocity, and referenced in symbol type to the CMD plot. 




Figure 5. The same as for Fig. 4 for a field further along the Stream 'D' structure (F37). Symbols are again highlighted which have 
photometric metallicities consistent with the [Fe/H] distribution of Stream 'D' in Ibata et al. (2007): —1.7 < [Fe/H] < —0.7. 
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Figure 6. The same as for Fig. 5 for stars in the Stream 'B' field, named 'M6' from Gilbert et al. (2006) and Koch et al. (2008). Symbols 
are again highlighted which have photometric metallicities consistent with the [Fe/H] distribution of Stream 'B' in Ibata et al. (2007). 
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Figure 7. The same as for Fig. 5 for stars in the Stream 'A' field, named 'M6' from Gilbert et al. (2006) and Koch et al. (2008). Symbols 
are again highlighted which have photometric metallicities consistent with the [Fe/H] distribution of Stream 'A' in Ibata et al. (2007). 



given structure are likely far more reliable than those mea- 
sured from the relatively low-S/N spectroscopy. 



3.2 Stream 'Cr' 

Stream 'C is the dominant stellar component at the po- 
sition of our spectroscopic masks F25/F26, exceeding the 
halo stars as well as stars likely to be foreground MW con- 



taminants (in the adopted velocity range of the M31 halo: 
<-150 km/s) by a factor of ~3x on average. In Fig. 2, the 
stars likely belonging to Stream 'C, and EC4 are highlighted 
(where EC4 member stars were preferentially inserted in the 
spectroscopic masks) allowing metallicity comparison to the 
fiducial globular cluster RGBs. To differentiate EC4 stars 
from the halo field, the velocities are also plotted against 
their radius from the EC4 center. Figure 2 shows that once 
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Figure 8. Velocity differences of stars lying in both fields F25 and F26. Stars identified to velocity regions likely associated to Stream 'Cr', 
Stream 'Cp', EC4, background M3f halo, and Milky Way foreground are identified. Stars are shown at their velocities from mask F25, 
with offset error-bars from mask F26. 



the cluster, EC4, stars are removed, a strong metal rich 
peak of stars at ~ —350 km/s dominates the stars kine- 
matically identified to exclude the Milky Way. The stars 
in this kinematic structure have average photometrically 
derived metallicity, [Fe/J?]=-0.74±0.19. As the metallicity 
of this kinematic substructure is very similar to that mea- 
sured for the total Stream 'C ([Fe/H]=-0.6) by Ibata et al. 
(2007), this is an excellent candidate for a kinematic de- 
tection of Stream 'C. Taking the clump of metal rich stars 
±2cr from the peak, we find nine stars with an average 
< v r >= -349.5ljg km/s, a Vr = 5.ljlf;| km/s, where 
the individual velocity errors are taken into account in a 
maximum-likelihood sense. This procedure is described ex- 
actly in Martin et al. (2007), although no iterative clipping 
is done in this case as the contamination levels from fore- 
ground Milky Way stars is much smaller. Briefly, using only 
the candidate Stream 'Cr' stars, a maximum-likelihood al- 
gorithm that explores a coarse grid of the (v r , a) space and 
searches for the couple of parameters that maximizes the 
ML function^ defined as: 

N 

ML(v r , cr) = ^log(— exp f - \ ^ )"} ) (1) 

i=i 

with N the number of stars in the sample, v r ,i the radial 
velocity measured for the i star, « err ,i the corresponding 
uncertainty and a tot = ^/cr 2 + v'j rr i . Using this definition 
of a allows to disentangle the intrinsic velocity dispersion 
and the contribution of the measurement uncertainties to 
these likelihood distributions^. These distribution functions 
are shown in Figure 9. 

1 There is an error in this expression from Martin et al. (2007) 
which is corrected here. 

2 An alternative way to determine v r and a is to use a = a' tot 
in equation (f) to measure the observed dispersion and then cor- 
rect from the mean velocity uncertainty, v err , such as cri t = 



The velocity is close to the systemic velocity of M31, - 
300 km/s, which might be expected if the stream were truly 
close to tangential as it appears to be in the imaging (most 
of its velocity being orthogonal to our measured heliocen- 
tric component). We cannot estimate a reliable mass for the 
progenitor from this single spectroscopic measurement, as 
the stars we have detected lie off-center from the Stream 'C 
peak, and as yet we have not measured the full extent of the 
stream. Nonetheless, we can at least place a constraint on 
the mass from the measured velocity dispersion in this field 
(§ 4). 

By taking the average halo profile of Ibata et al. (2007) 
at this projected radius, ~9 true halo stars are expected in 
fields F25/F26 (all stream structures removed from consid- 
eration), assuming all possible candidate RGBs have been 
observed in the two overlapping DEIMOS pointings, which 
they have. We find 26 halo stars are detected in the velocity 
region Vhel < —150 km/s excluding the MW, and after re- 
moving the high confidence EC4 stars (lying within 2 core 
radii) from consideration as they were added selectively to 
the mask in addition to the randomly selected halo stars. 
If 9 stars are associated to Stream 'Cr' and as we shall see 
another 7 stars are associated to Stream 'Cp', there are 10 
candidate halo stars found in the sample, in good agree- 
ment with the average prediction. Integrating the windowed 
a v ~ 125 km/s Gaussian of the halo, we find 3% chance 
that a star lies in one of the three 10 km/s velocity bins 
encompassing this kinematic structure. At least one, but 
unlikely more than two of the candidate Stream 'Cr' stars 
will be unrelated halo. Therefore neither the metallicity or 
velocity dispersion is likely to be heavily biased by unrelated 
halo stars. We note that << 1% Galactic contamination is 
expected at ~ —350 km/s, from our own characterization 



\J a 1 + Verr,i 2 ■ Parameters obtained in this way are similar to 
those given in the text. 





Figure 9. The kinematics of Stream 'C Upper panels: The relative likelihood distribution (taking into account the measurement 
errors) of < v r > and u Vr for velocities of Stream 'Cr' stars when marginalizing with respect to the other parameter. The thin dashed 
lines correspond, from top to bottom, to the parameter range that contains 68.3%, 95.4% and 99.73% of the probability distribution 
(1, 2 and 3x <r uncertainties), revealing that Stream 'Cr' has a none-zero velocity dispersion at the lex level, and a remarkably small 
<j Vr . The peak values and la range are < v r >= — 349.5^'g and a Vr = 5.1^2 \- Lower panels: The same for Stream 'Cp' (with 6 
candidate members), showing a none-zero velocity dispersion at the > 3(T level. The peak values and la range are < v r >= —285. 6^2 
and a Vr = 4.3^'^. Note that both streams have maximum likelihood distributions close to symmetric about the solution in projection. 



of the MW population in our spectroscopic fields, from the 
Gilbert et al. (2006) analysis of MW dwarfs in their M31 
spectroscopy, and from the Besangon Galactic populations 
model (as described in Ibata et al. 2005, 2007 and Chapman 
et al. 2006). 

3.3 Stream 'Cp' 

In analyzing the velocity distribution of halo stars in field 
F25/26, we move on from the strong kinematic peak of stars 
at -349 km/s which we have identified with Stream 'Cr', and 
notice in figure 2 a kinematic association of stars at Vhci ~- 
286 km/s showing a remarkably small dispersion. Figure 2 
also plots the radial distance of the stars in the field from the 
cluster EC4, where it is apparent that this spike merges with 
likely EC4 member stars. There are five stars lying within 
one core radius of EC4 (<30pc), two borderline members 
between 2-3 core radii, and 6 candidate stream stars at such 
large distance that they are unlikely to be directly associated 
to EC4. 

We note that one of these stars is likely either a blue 
supergiant star at M31 distance or is a Milky Way contami- 
nant, although the equivalent width of the Nal doublet is re- 
markably small for a Milky Way dwarf (e.g., Guhathakurta 
et al. 2006, Martin et al. 2007), with the caveat that this 
example is very blue in colour where it is not clear Nal is a 
good discriminant (e.g., Koch et al. 2008). Since no photo- 



metric metallicity can be derived for this star, we will not 
consider it further in our analysis (the star is Stream 'Cp', 
00:58:24.32+38:04:29.9). 

A tight range in [Fe/H]=-1.26±0.16 is observed in the 
six unambiguous Stream 'Cp' stars, very close to what we de- 
rive from the Ibata et al. (2007) foreground subtracted pho- 
tometry dataset for the offset Stream 'Cp' region, [Fe/H]=- 
1.1. We propose that in addition to Stream 'Cr', we have 
also kinematically detected this lower contrast Stream 'Cp' 
in our fields F25/F26. It is a matter of debate how the 
two borderline EC4 stars are treated; starting from the 
Stream 'Cp' standpoint we find no valid reason to reject 
them from the stream sample, as they lie well within the 
velocity window defined by the stars at much larger ra- 
dius (see Fig. 2), although we quote our results with and 
without them. Again, the systemic velocity close to that 
of M31 (vh e ;=-285.6±1.2 km/s with 6 members, Vh e i~- 
285.0^ 8 kni/s with 5 members) would be consistent with 
the expected properties of a tangential stream. Using the 
maximum likelihood technique we calculate a true veloc- 
ity dispersion of o\, iCorr = 4.3^^4 km/s with 6 members 
(<J v ,corr = 5.li^5 W s with 5 members). The likelihood 
distribution functions are shown in Fig. 9. 

The stacked spectra of Stream 'Cr' and Stream 'Cp' are 
shown in Fig. 10, emphasizing the clear difference in spectro- 
scopically derived [Fe/H] (-0.8 versus -1.4) between the two 
kinematic peaks, in good agreement with the photometric 
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Figure 10. The stacked spectra of Stream 'Cr' and Stream 'Cp', weighting by the inverse variance in the continua, emphasizing the clear 
spectroscopic difference in [Fe/H] (-0.8 versus -1.4) between the two kinematic peaks, in good agreement with the photometric [Fe/H]. 



[Fe/H] quoted in Table 1 (despite the fact that individual 
spectra are low signal-to-noise and show a large spread in 
[Fe/H]). With confidence that we have truly detected two 
different, superposed stream components through their off- 
set kinematics and metallicities, we revisit the Stream 'C 
region from fbata et al. (2007). In Fig. 11 we show a zoomed- 
in region around Stream 'C, divided in [Fe/H] into two non- 
overlapping ranges (0.0 to -0.7, and -0.7 to -1.7). This divi- 
sion makes it clear that two physically separate structures 
are present, and we can attempt to separate the luminosities 
of the two components (see § 4). 

With conservatively only 5 member stars (the 6th not 
obviously being an RGB star), and a metallicity range well 
within that expected for M31's overall stellar halo ([Fe/H~ 
-1.4 ± 0.2 from Chapman et al. 2006; Fe/H~-1.2 along the 
minor axis in Kalirai et al. 2006 or Fe/H~-1.5 over the same 
minor axis fields in Koch et al. 2008), we should first consider 
how likely this kinematic peak is to be distinct from the 
smooth halo component. It is also of interest to demonstrate 
that these stars are not likely to be far-flung members of 
EC4. 

These five stars at large radii from EC4 (0.7 to 10.3 ar- 
cmin, or 150 pc to 2220 pc), chosen by chance in the spectro- 
scopic masks are unlikely to be bound members of EC4 itself. 
Firstly, these stars would represent 5-74 core radii of EC4 
(30 pc). It is implausible that stars this far from the center 
belong to the cluster unless it is strongly disrupted, or is in 
fact only the cold core component of a more diffuse dwarf 
galaxy with an outer second component of stars. EC4 does 
not show any obvious signs of disruption, but given the faint- 
ne ss of EC4, it i s diffic ult to tell from the HST/ACS image 
in lMackev et al.l (|2006l ) whether there are subtle signatures 
of disruption. It can also be seen in Fig 2 that the metallici- 
ties of these stars appear to be marginally richer than EC4, 



although photometric errors from EC4 star crowding could 
easily account for these differences. 



It is more difficult to differentiate these stars from M31 
halo stars than it is for Stream 'Cr' stars, since the peak of 
the halo velocity distribution lies at ~-300 km/s, although 
the very broa d velocity dispersion ^ 125 km/s at this pro- 
jected radius (|Chapman et al.ll2006l ) makes it less likely to 
find such a strong spike of stars at -286 km/s. We estimate 
the chance association specifically as follows. We assume we 
know nothing about EC4 and that we have an ensemble of 
stars which are members of the M31 halo. For this purpose 
we assume conservatively that a halo star is any star con- 
sistent with the halo CMD, and to coarsely remove Galactic 
contaminants, lying within —550 km/s < v r < —150 km/s, 
or roughly ±2cr with a window clip appropriate to the 
Galactic contaminant distribution in Fig. 2. Removing the 
stars specifically targeted to lie in EC4, there are 26 stars 
which satisfy these criteria. There is a 3.2% chance that 
any given halo star will lie in the 10 km/s window centered 
on the five unambiguous Stream 'Cp' stars systemic veloc- 
ity (-285 km/s). This is a conservative assessment since the 
chance would be much lower for a window offset from the 
~-300 km/s peak of the M31 halo distribution. However, the 
chance that 5 stars out of 26 lie in this window is extremely 
small (irrespective of how the window is defined, the prob- 
ability is consistently below < 10 4 ), given the broad halo 
<j v . This estimate is even more conservative as there is in 
fact a substantial contribution from the Stream 'Cr' stars 
in this halo field - a better estimate of the total underly- 
ing halo stars would be ~ 16. We can conclude that these 
5 stars represent a rare kinematic spike in a smooth halo 
distribution. 
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Figure 11. The Stream 'C region shown zoomed-in (from Fig. 20 of Ibata et al. 2007, which details the definitions of cutout regions etc.) 
with slices in [Fe/H] (-0.7 < [Fe/H] < 0.0 - left panel; -1.7 < [Fe/H] < -0.7 - right panel) allowing the spatially offset Stream 'Cr' 
(left panel) and Stream 'Cp' (right panel) components to clearly be seen. Our kinematic separation of these two components appears 
to be reflected in physically distinct (but somewhat overlapping) regions. The comparison also highlight the much more metal poor 
Stream 'D' in the right panel. 



3.4 The continuations of Stream 'C 

As described in § 2, upon detecting cold kinematic peaks 
plausibly attributed to Stream 'C, we observed an addi- 
tional Keck/DEIMOS field (F36) further North along this 
structure. The properties of this field are shown in Figure 3. 
While no prominent peaks are found in the velocity distri- 
bution of this northern field, a concentration of five metal- 
rich ([Fe/H]~-0.7) stars is observed at an average velocity 
of -350 km/s, possibly attributable to the same Stream 'Cr' 
structure, and thereby showing no obvious velocity gradient. 

In the same field (F36) we also search for the more 
metal-poor Stream 'Cp'. We find a clump of four similarly 
metal-poor stars, <[Fe/H]>=-1.2, at an average velocity of 
-246 km/s, offset by ~ +40 km/s from the Stream 'Cp' and 
EC4 peak in fields F25/26. 

Orbit models of these streams will be presented in a 
future paper, while Fardal et al. (2008) discuss how these 
stream-like structures could conceivably be related to the 
Giant Southern Stream. 

3.5 The extended star cluster, EC4 

Figure 2 shows that stars targeted in the cluster, EC4, 
have clearly been kinematically identified with a distribu- 
tion of velocities centered at Vhei=-285 km/s. A detailed 
study of EC4 (Collins et al. in preparation) suggests a small 
(~ 3 km/s) resolved velocity dispersion, marginally consis- 
tent with a dark matter dominated system. The EC4 stars 
have [Fe/H] p hot=-1.4±0.1, although a spectroscopic esti- 
mate of the [Fe/H]=-1.6±0.15 0. As these stars were pre- 
selected to lie in EC4, we assume that at least the 5 stars 

3 Mackey et al. (2006) find [Fe/H] pfeot =-1.84 for EC4 from HST 
photometry and simultaneous fitting of the RGB and Horizontal 
Branch. This difference is explored in Collins et al. (in prepara- 
tion). 



lying within 1 core radii of the EC4 center (and likely the 
two stars at 2-3 core radii) can be removed from the sur- 
rounding M31 halo sample for our statistical analysis of the 
Stream 'C in previous sections. We have clearly identified 
the systemic velocity of EC4 as being compatible with the 
Stream 'Cp' kinematics. 

3.6 Stream 'D' 

Based on the example of Stream 'Cr' and Stream 'Cp', we 
are motivated to search for a narrow velocity peak in the 
case of Stream 'D'. However Fig. 4 showing the CMD, radial 
velocities and metallicities (as for the previous streams) does 
not reveal any obvious detection of stars in this stream. We 
proceed by comparing the expected metallicity from Ibata 
et al. (2007), -1.7 < [Fe/H] < -0.7, with any stars in 
the "halo" sample (culled from the velocity and EW(Nal) 
cuts) which could be the Stream 'D'. We highlight all stars 
in the CMD which could conservatively be consistent with 
the Stream 'D' median [Fe/H]=-1.2. There are no obvious 
kinematic spikes within these colour-selected stars, as shown 
in Fig. 4 (right panel). However, there is an isolated group 
of two stars within a 10 km/s bin at -405 km/s. Two stars 
at ~-400 km/s are somewhat unexpected (8% chance) given 
the halo velocity dispersion at a projected radius of 35 kpc 
(Chapman et al. 2006). There are no better candidates for 
the Stream 'D' than this pair of RGB stars, but we cannot 
confidently separate Stream 'D' stars from spheroidal halo 
stars. 

In the DEIMOS field placed further along Stream 'D' 
(field F37), an kinematic association of 5 stars stands out 
again at ~400 km/s (as with field F7 above), however only 3 
have inferred [Fe/H] within a range consistent with the pho- 
tometric properties/ If these stars represent the Stream 'D', 
there is also no measurable velocity gradient detected, as 
found for Stream 'Cr'. In Table 1 we present the velocity 
dispersion measured from the 5 plausible Stream 'D' stars 
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combined from both fields (a Vr =4.2 km/s). However, the 
low contrast of Stream 'D' relative to the background M31 
halo, together with our inability to distinguish with con- 
fidence a kinematic identification mean that these results 
are tentative. Much larger numbers of spectroscopic mea- 
surements along Stream 'D' are required in order to reliably 
detect a coherent structure in velocity. 



3.7 Streams 'A' & 'B' 

The other two streams lying perpendicular to the minor axis 
presented in Ibata et al. (2007), Stream 'A' and Stream 'B' 
lie at 120 kpc and 80 kpc respectively. Both of these streams 
have serendipitous spectroscopic pointings lying in their 
edge regions from Gilbert et al. (2006) and Koch et al. 
(2008), named fields M8 and M6 respectively in their nomen- 
clature (see Ibata et al. 2007 for placements of these spectro- 
scopic pointings in the wider M31 halo map). We present the 
CMDs, velocity histograms and metallicities for Stream 'B' 
and Stream 'A' here for analysis (Figs. 6 & 7). 

In field M8 (Stream 'A') there are only 4 stars with 
velocity measurements attributable to the M31 halo. The 
clump of three stars at —172 km/s have an [Fe/H]=-1.3 on 
average, very similar to the statistical measurement of the 
[Fe/H]~ -1.3 in Stream 'A' (Ibata et al. 2007). We further 
suggest, along Ibata et al. (2007) and Koch et al. (2008), that 
these 3 kinematically identified halo stars are more likely to 
be associated to the stellar halo of M33 (the Triangulum 
galaxy), given the M31 halo velocity dispersion (Chapman 
et al. 2006) shown in Figs. 6 & 7. It is therefore worth con- 
sidering that this Stream 'A' structure might actually be dis- 
rupted remains of a satellite in M33's halo. The dispersion 
of these three stars, 14.5 km/s, cannot easily be deconvolved 
for measurement errors using the maximum likelihood ap- 
proach. Instead, since the errors are similar for all three 
stars, we write as in footnote(l) a VtCorr = — o\ natT ) 

= VTl4.5 2 - 7.2 2 ) = 12.5 km/s. 

We carry out the same procedure as with the other 
streams, identifying candidate stars in the CMD which 
are consistent with the average metallicity found in Ibata 
et al. (2007) for the streams in question. In field M6 
(Stream 'B'), a kinematic peak of stars lying at ~-330 km/s 
with <[Fe/H]>~ —1.0 represents a reasonable candidate for 
this stream. Notably, this metallicity is very close to that 
estimated for Stream 'B' in Ibata et al. (2007), and further, 
both the metallicity and the velocity distribution in this field 
depart significantly from the average found in Koch et al. 
(2008) for the outer halo. Again following the logic of our 
discussions in the streams 'C and 'D', the RGB overden- 
sity in this field attributed to the stream should statistically 
result in the bulk of stars kinematically detected lying in 
the stream structure. Removing these five relatively metal- 
rich stars from the halo sample leaves six more metal-poor 
stars which would represent the surrounding M31 halo at 
this radius. As with the Stream 'A' case above, we simply 
estimate an intrinsic dispersion directly from these five can- 
didate Stream 'B' stars, using their average measurement 
error, of a Vr = 6.9km/s. 



4 THE STREAMS IN THE CONTEXT OF THE 
DWARF SPHEROIDAL M/L TO L 
RELATION 

It is of interest to ask how the properties of these streams 
would compare to other M31 satellites if treated as dwarf 
remnants. As described in Ibata et al. (2007) the light in 
polygonal regions surrounding each stream was integrated, 
with background corrections applied. Here we go further 
to emphasize the metal-poor region of Stream 'C defined 
as the V-band light with -2.0 <[Fc/H]< -1.0, and fur- 
ther assuming (through the luminosity ratios of the offset 
portions of the metal-poor and metal-rich components, de- 
scribed above) that 1/10 of this metal-poor light belongs 
to Stream 'Cp' while the other 9/10 belongs to Stream 'Cr'. 
The resulting luminosities are summarized in Table 1. There 
is of course an uncertainty in the Mv estimates, since the 
streams are terminated by the edge of our M31 halo imaging 
to the North. We will therefore simply assume here that the 
My estimates are lower limits. 

The width of the streams can provide constraints on the 
mass of the progenitor, since debris from more massive satel- 
lites is likely to produce wider debris streams that spreads 
more rapidly along the orbit with time. We measure the 
widths of all our streams by taking the minor axis profile in- 
tegrated over the full extent where the streams are detected. 
The streams are not always exactly orthogonal to the mi- 
nor axis, however they are close enough (and in any case 
not always well defined) that any broadening introduced by 
this simple procedure should be minimal. We then fit Gaus- 
sians to the profile, subtracting the local background from 
regions on either side of each stream. We quote the FWHM 
in Table 1. 

Johnston et al. (2001) present simple analytic scalings 
for the width and length of debris streams, with the main 
assumption that the progenitor is supported by random mo- 
tions. In this case, the measured fraction, s = w/R, of the 
width to the radial distance to the stream is related to the 
mass m of the satellite through the relation 

s = (Gm/[v 2 circ R pm ]) 1,A . 

Font et al. (2006) estimate the progenitor mass of the 
Giant Southern Stream (GSS) in M31 using these mea- 
surements (w, R) and an orbit model suggesting R per i = 
3 - 4.5 kpc, finding 1.0 - 1.6 x 10 8 M Q . With an updated 
luminosity for the GSS from the wider survey of Ibata et 
al. (2007), 1.5 x 10 8 M Q , the M/L ~ 1 M Q /L . We note 
however that if the models of the GSS progenitor as rota- 
tionally supported (Fardal et al. 2008; Mori & Rich 2008) 
are correct, this would invalidate the Font et al. estimate of 
progenitor mass. 

We do not yet have sufficient information to model the 
orbits of these five fainter streams in M31, and thus our esti- 
mates of the progenitor masses are slightly less secure than 
the GSS estimate in Font et al. (2006), with a linear de- 
pendence on the uncertain R per i- We place limits on R per i 
by assuming that we are seeing these streams at apocen- 
tre, and that orbits in cosmological simulations have aver- 
age R P eri/Ra P o ~ 0.2 0.25 (van den Bosch et al. 1999; 

Ghigna et al. 1998; Benson 2005). These limits are listed in 
Table 1 along with the other model parameters (and taking 
for M31, u clrc =260 km/s). 
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Figure 12. Comparison of the mass-to-light ratios (M/L) and light, My for the streams in M31, the Sagittarius stream in the Milky 
Way (Majewski et al. 2003), and compared to the faint dwarf galaxies from both the MW (Simon & Geha 2007; Martin ct al. 2007) and 
M31 with central velocity dispersion estimates. The procedure for measuring the M/L is described in the text. Stream M/L values are 
shown as bars connecting the structural mass estimate to the a v estimate. A factor of two uncertainty is shown for the a v mass since we 
are likely observing the streams at especially cold points between turning, and the true progenitor mass is likely 2-4 times larger. And XI, 
And XII, and And XIII (Collins et al. in preparation) are all shown as upper limits, since their velocity dispersions are all unresolved by 
their measurements — the Icr likelihood contour is used to set a tentative limit. AndXV & AndXVI measurements come from Letarte et 
al. (2008). The solid lines are curves of constant dark matter halo mass (1, 2, 4, 8x 10 7 Mq from bottom to top), assuming a stellar 
mass-to-light ratio of 2.5 Mq/Lq. 



The same procedure can also be applied to the Sagit- 
tarius stream (Sgr), taking parameters from Majewski et 
al. (2003), with R per i = 12 kpc, wfwhm=4 kpc, v circ = 
220 kpc, we calculate a mass of 5.5 x 10 9 M©. Majewski et al. 
(2003) derive a mass from kinematics of Sgr of 5.8xl0 8 Mq 
with a M/L=25 (Mv = -13.27). 

The velocity dispersions of the streams can also be used 
to constrain the progenitor masses, although perhaps with 
even less accuracy. Generally, the dispersion in stream de- 
bris should decrease over time (Helmi & White 1999b). If 
these streams were very young, we wouldn't expect dynam- 
ical cooling to be signicant yet. Tidal interactions with dark 
matter substructure in the halo may also not have had suffi- 
cient time to significantly heat the streams if they are rela- 
tively young (Ibata et al. 2002, Johnston et al. 2002). Helmi 
& White (1999b) suggest the velocity dispersion should vary 
most signicantly in an oscillatory manner as a function of 
radial orbital phase. 

Font et al. (2006) also estimate the progenitor mass of 
the GSS from a single cr Vr — 15 km/s lying between apoc- 
entre and pericentre along the stream, with a lower limit of 
10 8 Mo, consistent with their estimate from structural/orbit 
properties. As predicted by Helmi & White (1999), the 
stream can become very cold in between the turning points 



with the velocity dispersion of the stream reaching values 
well below the central dispersion of the satellite, and as small 
as a /(Jo ~ 0.5. 

In a similar manner, we can estimate the M/L of the 
streams presented in this contribution. For all the streams, 
we can apply the methodology of Font et al. (2006) directly 
as we have plausible kinematic detections in each case, along 
with reasonable constraints on the stream widths and mor- 
phologies. We expect the velocity dispersion of the progeni- 
tor to be as large or larger than the intrinsic value estimated 
for the stream and assume mass follows light (Richstone & 
Tremaine 1986). For Stream 'D', we have taken the combined 
a Vr from both spectroscopic pointings (Table 1), assuming 
that the ~-400 km/s stars are the most likely members of 
the structure. For Stream 'A', we assume the three halo stars 
at —172 km/s represent the stream (although as discussed, 
it is ambiguous whether it is a structure associated to M31 
or M33). The results of our mass estimates are presented in 
Table 1, and plotted on Fig. 12, where we highlight a fac- 
tor of two uncertainty given that we are likely observing the 
streams at especially cold points between turning. 

However, we caution that these M/L ratios have sig- 
nificant uncertainties attached. The progenitor may not be 
completely disrupted, the velocity dispersion could be a very 
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poor mass estimate (more typically a lower limit) depending 
on where in the orbit and evolution the stream is, and from 
the truncation of the image to the North, we already know 
the observed luminosity may not be representative for the 
entire streams. 

The comparison of the streams with the Milky Way 
satellites and all published M31 satellites is presented in 
Figure 12. Stream M/L values are shown as bars connect- 
ing the structural mass estimate to the a v estimate. The 
Mateo (1998) relation between M/L and the luminosity of 
a dwarf galaxies in the Local Group has also been plot- 
ted on Figure 12. The relation can be understood physi- 
cally as more massive dwarfs retaining more of their gas 
(and therefore arriving at z — with a lower M/L ratio), 
while lower mass halos more easily expel their gas and form 
smaller numbers of stars. All of our newly constrained M31 
streams are in agreement with the relation within reasonable 
errors, which is somewhat surprising given our limited abil- 
ity to constrain the masses of the streams. Interestingly, the 
Sagittarius stream lies significantly above the relation over 
all plausible mass estimates. Many of the new M31 dSphs, 
along with most MW dSphs fall on the relation, although 
some clearly appear to be outliers, as are many of the new 
faint Milky Way dSphs (Martin et al. 2007; Simon & Geha 
2007). This suggests that at the low mass end, a range of pro- 
cesses beyond simply the feedback in winds which explains 
galaxies as massive as our stream progenitors may truncate 
star formation (e.g., Ricotti & Gnedin 2005), or else the nu- 
merous model assumptions are failing when reaching these 
faint limits. It will be of interest to obtain improved con- 
straints on these M31 streams and model their orbits, to see 
if progenitor mass estimates continue to keep them on the 
Mateo relation. 



5 DISCUSSION 

Our spectroscopic survey of the M31 faint streams has 
yielded an encouraging initial census of the kinematics, how- 
ever it is clear from our results that significant efforts with 
a 10 metre telescope are required to usefully constrain the 
kinematic properties of the streams for modeling. At least 
in the case of the spatially overlapping Stream 'Cr' and 
Stream 'Cp', it is remarkable that we have been able to 
clearly distinguish these structures by kinematics. Indeed, it 
was not entirely clear from initial inspection of the imaging 
in Ibata et al. (2007) that there were two different streams 
in this vicinity at all. Only careful inspection of the imaging 
divided in slices of metallicity reveals two structures slightly 
offset spatially. 

The cluster, EC4, lies in a region where the metal-poor 
Stream 'Cp' has roughly 50% the stellar density of the metal- 
rich Stream 'Cr', although our spectroscopy reveals that EC4 
is likely related to Stream 'Cp' with Stream 'Cr' overlapping 
only in projection. Could Stream 'Cp' actually be the debris 
from disrupted EC4 material? The integrated luminosity of 
Stream 'Cp' within the MegaCam survey is comparable to 
a small dwarf galaxy like And XV or And XVI (Letarte et 
al. 2008), My ~ —9.5, which would suggest the baryonic 
matter mass loss of EC4 (My = -6.6, Mackey et al. 2006) 
would be dramatically larger than its current intact mass. 
No distortion of the EC4 isophotes is found in the HST 



imaging of Mackey et al. (2006) (Tanvir et al. in prep), al- 
though the faintness of EC4 means this is unlikely to be a 
good test of ongoing mass loss or tidal distortion. We have 
however noted that the outer stars in EC4 show a statisti- 
cally significant velocity shift from the inner EC4 stars, sim- 
ilar to the Stream 'Cp' stars in the field surrounding EC4. 
This could happen for instance if EC4 were disrupting in a 
stream along the line of sight. Regardless, it is likely that 
Stream 'Cp' and EC4 are at least related by their kinematics 
and metallicities, EC4 possibly representing an intact sys- 
tem carried along in the disrupted progenitor represented 
by Stream 'Cp'. If EC4 is dark matter dominated, we have 
in fact detected the very first sub-sub-halo (i.e. a galaxy 
that was bound to a satellite galaxy), possibly explaining 
its small (r c =30pc) size. From a ACDM cosmological point 
of view, the LMC and SMC should also have such sub-sub- 
halos with L ~ 10 7 -10 8 L , but we find none (and here 
we cannot invoke tidal disruption of these systems, because 
they seem to be falling in for the first time) . Of course it also 
remains the possibility that EC4 has nothing to do with the 
Stream 'C structure at all. Without precise distance infor- 
mation, it is difficult to rule this out completely. 

We also ask whether these stream structures could be 
related to the giant southern stream imaged in Ibata et 
al. (2001, 2007)? As we have noted, the metallicities of 
various streams all differ, only Stream 'Cr' and Stream 'B' 
being even close to the metallicity of the core region 
of the giant southern stream (Ibata et al. 2001, 2004, 
2007; Guhathakurta et al. 2006). However, in the CFHT- 
MegaCam survey, the outer region of the GSS are more 
metal-poor and have a metallicity that is similar to that 
of some of the more metal-poor streams. Nonetheless, the 
general impression from the low dispersion of the streams in 
kinematics, their physical thickness, and varying (but nar- 
row) metallicities makes it hard to reconcile with a scenario 
whereby these stream structures as an ensemble are related 
to the giant southern stream. 

Fardal (2008, and private communication) has modeled 
the stream resulting from a GSS progenitor that is flattened 
and rotating like a disk, building on the constraints and 
models from Fardal et al. (2007) and Gilbert et al. (2007). 
Because it is on such a radial orbit, when it reaches pericen- 
ter very close to M31, part of the progenitor could be on the 
opposite side of M31 to the rest, depending on the orienta- 
tion of the disk. This means that it starts orbiting M31 in 
the opposite direction to the rest and leaves debris in differ- 
ent physical locations than the main stream. The disk-like 
kinematics results in caustic structures that appear similar 
to streams or arcs. In this model, the new streams from Ibata 
et al. (2007), and herein, would be shells from this counter- 
orbiting part of the stream. Assuming a large metallicity 
gradient in the progenitor, they would have a much lower 
metallicity than the main stream since they come from an 
outer part of the progenitor. Whether or not this specific 
model is correct, the general idea that the progenitor was 
physically quite large, and passed extremely near the center 
of M31, means that debris could get thrown out in all direc- 
tions. Given the metallicity structure of the progenitor, it is 
plausible that some of this debris could have distinct metal- 
licities but be ultimately related to the same progenitor. 
While somewhat implausible for the reasons stated above, 
it remains to be seen if the specific kinematics and metal- 
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licities of our new observations can be reproduced in such a 
model. 

Whereas the photometric profiles could only remove 
Galactic contamination and stars belonging to stream sub- 
structures statistically, we can explicitly remove stars be- 
longing to the streams (and the Milky Way) by their kine- 
matics and assess the underlying M31 stellar halo density 
from 30-120 kpc on the minor axis. While our resulting mea- 
surement has so few stars as to be highly uncertain statisti- 
cally, it does reveal the general power of kinematic analysis 
of the M31 halo population. 

It is remarkable the extent to which these kinematic 
substructures projected on the minor axis dominate the halo 
star statistics in these fields. In the fields studied, they rep- 
resent ~ 2/3 of the candidate halo stars, revealing that the 
photometric minor axis profile from Irwin et al. (2005) and 
Ibata et al. (2007) is significantly flattened by such struc- 
tures. We are led to the likely conclusion that stellar halos 
are made up of multiple kinematically cold streams, perhaps 
even to the extent proposed by Bullock & Johnston (2005) 
(see also Bell et al. 2008). 



([Fe/H]=-1.4) which suggest it is related to the more metal- 
poor stream Stream 'Cp'. EC4 could be the progenitor of the 
metal-poor Stream 'Cp' (somewhat unlikely given the appar- 
ent stellar mass difference between the stream and EC4), or 
it may simply be a structure carried along by the disrupted 
stream progenitor. In this case, and if EC4 has a sizable 
dark matter component, we have in fact detected the very 
first sub-sub-halo (i.e. a galaxy that was bound to a satellite 
galaxy), possibly explaining its small (r c =30 pc) size. 

• By explicitly removing stars belonging to the streams 
by their kinematics we can assess the underlying M31 stellar 
halo density and metallicity on the minor axis. This con- 
trasts the purely photometric approach where Galactic con- 
tamination and stars belonging to stream substructures can 
only be removed statistically. Our resulting halo measure- 
ment has so few stars as to be highly uncertain statistically, 
however, it does reveal the general power of kinematic analy- 
sis of the halo population for future endeavors. The fraction 
of background halo stars in these stream fields suggests the 
conclusion that stellar halos are largely made up of multiple 
kinematically cold streams. 



6 CONCLUSIONS 

In conclusion, we have conducted a Keck/DEIMOS spec- 
troscopic survey of five stellar streams, recently uncovered 
through deep imaging observations of the halo. 

• We have uncovered a kinematic substructure at Vh e i=- 
349.5±1.8 km/s from a spectroscopic field lying in the Ibata 
et al. (2007) Stream 'C. The cold component has a Vr = 
5.1 ± 2.5 and a narrow range in [Fe/H]=-0.7±0.2, which we 
propose represents a metal-rich component, Stream 'Cr'. 

• We have uncovered a second kinematic substructure 
in the same field as Stream 'Cr' at Vh c i=-285.6±1.2 km/s 
with o~ Vr — 4.3l}'4 km/s (non-zero at >3<r confidence in- 
terval) and a narrow range in [Fe/H]=-1.3±0.2, which we 
propose represents a metal-poor stream, Stream 'Cp'. We 
demonstrated that this kinematic Stream 'Cp' has a coun- 
terpart in a spatially offset metal-poor region of Stream 'C 
in Ibata et al. (2007). 

• We plausibly detect both Stream 'Cr' and Stream 'Cp' 
at a position ~30kpc further north along the structure, with 
no detectable velocity gradiant for Stream 'Cr', and a mea- 
sured velocity gradient of ~40 km/s for Stream 'Cp'. 

• We were unable to identify kinematic substructure un- 
ambiguously associated to Stream 'D' from our serendipitous 
spectroscopic pointing, however subsequent spectroscopy 
well centered in the Stream 'D' identifies a likely cold kine- 
matic structure which has a viable counterpart in the 
serendipitous pointing. We propose a kinematic detection 
of Stream 'D' at Vh c i=-390.5 km/s with a Vr = 4.2 km/s. 

• Spectroscopy near the edges of Stream 'A' and 
Stream 'B' suggest a likely kinematic detection for 
Stream 'B' with Uhei ~ —330 km/s, a v , corr ~ 6.9 km/s, and 
a kinematic detection of Stream 'A' at v^d ~ —172 km/s, 
cr Vr ~ 12.5 km/s. Neither spectroscopic pointing in these 
streams is ideally placed, and additional spectroscopic ob- 
servations are well motivated to further constrain the kine- 
matics of these structures. 

• The extended cluster EC4 lies in the Stream 'C 
region, with kinematics (vhoi=-285 km/s) and metallicity 
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Table 1. Properties of DEIMOS fields in the five M31 stellar streams 'A','B','Cp','Cr','D'. 



field 


a (J2000), 8 (J2000) 


IV, stream 


a 

° stream 


< [Fe/H] > b 


^stream 




M stream 


^stream 








(km/s) 


(km/s) 


(kpc) 


(kpc) 


xlO 7 M 


xlO 6 L Q 


stream A (M8) 


01:14:01.37 +32:31:00.9 


-172.2 


12.5 


-1.3 + 0.3 (-1.3) 


7.5 


24 


10/7 


2.3 


stream B (M6) 


01:32:14.64 +33:12:25.4 


-330.1 


6.9 


-0.8 + 0.2 (-0.6) 


5.0 


16 


7/2 


10.0 


stream Cr (F25/F26) 


00:58:22.02 +38:04:05.9 


-349.5 


5.1 ± 2.5 


-0.7 + 0.2 (-0.6) 


6.8 


12 


30/1 


12.6 


stream Cp (F25/F26) 


00:58:22.02 +38:04:05.9 


-287.3 


4 3+ 1 ' 7 


-1.3 + 0.2 (-1.1) 


8.5 


11 


70/1 


1.4 


stream Cr pos2 (F36) 


01:00:38.00 +38:45:37.0 


-350 


n/a 


-0.7 + 0.2 (-0.6) 


n/a 








stream Cp pos2 (F36) 


01:00:38.00 +38:45:37.0 


-246 


n/a 


-1.2 + 0.2 (-1.1) 


n/a 








stream D (F7) 


00:54:55.02 +39:43:55.3 


-390.5 


4.2 


-1.1 + 0.3 (-1.2) 


8.2 


6 


213/1 


9.5 


stream D pos2 (F37) 


00:57:34.00 +39:49:12.0 


-390.5 


4.2 


-1.1 + 0.3 (-1.2) 


n/a 








EC4 (F25/F26) 


00:58:15.50 +38:03:01.1 


-282.4 


~ 10 


-1.4+0.1/-1.69 


n/a 









a Velocity dispersions, estimated through a maximum likelihood analysis taking into account the measurement errors in the velocities 

(or in the case of Stream 'A', Stream 'B', Stream 'D', a subtraction in quadrature of the measurement error). 
6 For certain streams (notably Stream 'A' and Stream 'D'), it is arguable that the [Fe/H] p hot measurements are estimated from sets of 
stars which may not actually be a kinematic detection of "the stream". For this reason we also quote the statistical [Fe/H] p h t estimate 
from Ibata et al. (2007) in all cases in brackets. In particular the Stream 'D' values are quoted for the combination of 6 stars in the two 

separated pointings along the stream, and therefore the numbers in the table are simply duplicated. 
c Stream widths derived from Gaussian fits (quoted FWHM) to the integrated profile in the region defined in Ibata ct al. (2007) (their 
Fig. 31). At the distance of M31 (785 kpc - McConnachie et al. 2005) 1 degree = 13 kpc. 
d Rperi estimated as 1/5 D 'stream as discussed in the text. 
e Mass of stream width estimated from structural parameters (first entry) and from cr v (second entry). 
f Lstream from Ibata et al. (2007), except for Stream 'Cp', Stream 'Cr', which are discussed in the text. 
9 For EC4, the spectroscopic estimate of the [Fc/H] = -1.6, lying midway between the estimates derived using CFHT photometry with 
Girardi isochroncs (-1.4) versus and the HST photometry with Dartmouth isochroncs (-1.8) as described in Collins et al. (in 

preparation) . 

Table 2. Properties of candidate M31 halo stars in stream Di region (field F7). 



a (J2000) 8 (J2000) » r (kms- 1 ) « err (kms" 1 ) [Fe/H] hot [Fc/H] V-mag I-mag 



00:54:18.02 


+39:43:32.4 


-449.7 


6.4 




-3.12 


20.26 


19.38 


00:54:26.77 


+39:45:58.2 


-150.3 


3.2 


-1.78 


-2.027 


22.31 


21 


00:54:31.71 


+39:44:54.1 


-346.7 


9.96 


-3 


-1.662 


22.36 


21.3 


00:54:35.93 


+39:43:59.2 


-402.7 


11.05 


-1.11 


-3.295 


22.84 


21.32 


00:54:39.48 


+39:45:05.4 


-169.4 


34.95 




-2.181 


21.7 


21 


00:54:40.93 


+39:43:08.8 


-268.4 


11.45 


-0.73 


-2.843 


23.7 


21.29 


00:54:47.03 


+39:44:25.9 


-370.8 


8.01 


-1.58 


-3.454 


22.72 


21.44 


00:54:50.30 


+39:43:36.5 


-237.3 


14.53 


-1.43 


-2.144 


22.61 


21.22 


00:54:53.48 


+39:43:11.3 


-274.6 


5.46 


-3 


-2.637 


22.66 


21.64 


00:54:57.75 


+39:43:32.5 


-195.9 


49.42 


-1.22 


-3.394 


22.79 


20.76 


00:55:00.87 


+39:43:37.5 


-254.2 


25.61 


-1.12 


-3.135 


23 


20.82 


00:55:10.40 


+39:46:58.7 


-162.1 


2.73 




-2.025 


21.46 


20.7 


00:55:12.94 


+39:43:50.3 


-151.5 


2.44 




-2.427 


20.43 


19.62 


00:55:13.27 


+39:43:19.4 


-211.5 


7.26 




-2.894 


23.69 


20.87 


00:55:14.67 


+39:44:52.5 


-217.8 


38.62 




-3.182 


21.55 


21 


00:55:22.81 


+39:45:25.5 


-414.3 


8.37 


-1.18 


-1.537 


22.86 


20.87 


00:55:27.67 


+39:43:25.5 


-244.2 


3.11 




-2.73 


21.42 


19.16 


00:55:34.22 


+39:42:51.8 


-212.3 


12.14 




-2.776 


23.52 


20.75 
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Table 3. Properties of candidate M31 halo stars in stream D2 region (field F37). 



a (J2000) 


5 (J2000) 


v r ( km s 1 ) 


V err ( km S 1 ) 


[Fc/H] phot 


[Fe/H] spec 


V-mag 


I-mag 


nn- fis 


1 QQ.ixn.nfi Q 
-hoy .ou.uo.o 


-ZUO.D 


Q QA 




1 3fizl 

-± .Oo4 


1 Q 7^ 

iy. 1 o 


1O.O0 


nn-^7-no fin 

UU.O f .U-Z.oU 


-hoy. 4/ .OD. u 


1 KC\ O 

-lou.y 


i o A c: 


1 1Q 

-i.iy 


n /ioi o 

U.4Z1Z 


nn 7 


m io 


nn- i o 


I QQ./I Q. /I O O 

-hoy .4y .4y .y 


-1 ( O.D 


ft /1 7 
0.4 / 


-1.0/ 


U.oOll 


nn nn 

^z.oy 


Z1.U0 


uu.o 1 .uy.yo 


-hoy .Oz.oz.O 


-ODD.y 


i c: 01 
id. yi 


1 o 

-i.y 


-l.ool 


O/l 1 1 
Z4. 1 1 


ZL . ID 


00:57:35.45 


+39:49:42.4 


-187.6 


24.21 


-1.53 


2.214 


23.24 


22.1 


00:57:37.07 


+39:47:50.2 


-153.3 


4.8 




-0.9179 


21.02 


19.73 


00:57:40.78 


+39:50:55.2 


-264 


32.52 


-0.54 


-0.5198 


23.74 


21.93 


00:57:41.37 


+39:48:33.1 


-197.6 


3.56 


-0.44 


1.305 


23.93 


22.07 


00:57:43.61 


+39:51:13.5 


-192.6 


8.33 


-1.18 


-0.3217 


22.77 


21.28 


00:57:47.86 


+39:50:35.0 


-400.6 


43.33 


-1.39 


-0.605 


23.28 


22.1 


00:57:48.24 


+39:48:36.6 


-390.5 


10.44 


-1.04 


-0.3352 


23 


21.01 


00:57:50.48 


+39:49:09.6 


-169.2 


4.03 


-0.95 


-1.017 


23.27 


21 


00:57:55.32 


+39:51:25.7 


-386.6 


6.08 


-0.74 


-0.2843 


23.53 


21.08 


00:57:56.65 


+39:50:48.5 


-385.6 


8.55 




-0.9452 


21.1 


20.1 


00:57:57.41 


+39:51:43.2 


-233.7 


55.17 


-0.84 


-1.161 


23.42 


21.86 


00:58:00.09 


+39:48:33.5 


-247.3 


7.99 


-0.84 


-0.3931 


23.33 


21.33 


00:58:06.18 


+39:51:15.1 


-258.6 


15.52 




-1.195 


21.44 


20.61 



Table 4. Properties of candidate M31 halo stars in stream Ci region (field F25/F26). 



a (J2000) 


6 (J2000) 


v r ( km s 1 ) 


V err ( km S 1 ) 


[Fe/Hl h . 

L / iphot 


[Fc/Hl 

L / J spec 


V-mag 


I-mag 


Deca a 


00:57:46.69 


+38:11:43.0 


-349 


6.3 


-1.22 


-1.28 


22.49 


20.47 


10.38 


00:57:48.22 


+38:11:48.6 


-343.6 


8.38 


-0.89 


-4.83 


23.29 


21.72 


10.3 


00:57:51.02 


+38:11:25.4 


-193.3 


5.69 


-0.55 


-1.75 


23.74 


21.83 


9.68 


00:57:52.27 


+38:08:34.6 


-289.3 


8.68 


-1.01 


-1.76 


23.24 


21.82 


7.19 


00:57:54.36 


+38:12:26.8 


-286.9 


5.21 


-1.29 


-0.63 


22.57 


20.98 


10.3 


00:57:56.30 


+38:06:18.1 


-356.5 


6.7 


-0.62 


2.75 


23.7 


21.61 


5.00 


00:58:05.23 


+38:08:41.4 


-223 


2.22 


-0.4 


-1.8 


24.27 


21.48 


6.02 


00:58:12.75 


+38:05:24.6 


-341.8 


14.37 


-0.7 


-0.97 


23.58 


21.51 


2.45 


00:58:13.99 


+38:06:18.1 


-295.1 


5.96 


-1.26 


-1.09 


22.99 


21.66 


3.30 


00:58:14.37 


+38:03:00.4 


-296 


5.29 


-1.48 


-1.84 


22.91 


21.65 


0.22 


00:58:14.74 


+38:03:00.8 


-277.3 


9.56 


-1.49 


-5.62 


22.77 


21.47 


0.15 


00:58:15.24 


+38:03:01.3 


-293.4 


11.55 


-1.39 


-1.89 


21.7 


20.99 


0.05 


00:58:15.32 


+38:03:09.7 


-371.5 


10.9 


-2.17 


0.63 


22.65 


21.57 


0.15 


00:58:15.47 


+38:02:58.9 


-267.4 


8.42 


-1.45 


-0.53 


22.41 


20.92 


0.04 


00:58:15.94 


+38:04:34.9 


-153.3 


40.28 


-2.45 


3.38 


22.78 


21.73 


1.57 


00:58:15.99 


+38:02:56.1 


-264 


7.74 


-1.21 


-0.69 


22.75 


21.3 


0.13 


00:58:16.00 


+38:02:22.5 


-280.3 


2.83 


-1.41 


-1.38 


22.06 


20.41 


0.65 


00:58:17.12 


+38:02:54.1 


-281 


2.9 


-1.36 


-0.65 


22.49 


20.86 


0.34 


00:58:17.16 


+38:02:49.6 


-290 


2.57 


-1.16 


-1.12 


22.68 


21.09 


0.38 


00:58:18.62 


+38:00:16.7 


-530.7 


9.81 




0.73 


24.18 


20.6 


2.81 


00:58:19.68 


+38:06:42.4 


-577.2 


11.45 


-0.77 


-2.2 


23.5 


21.11 


3.78 


00:58:20.25 


+37:59:46.6 


-183.4 


19.86 


-0.56 


-0.96 


23.84 


21.65 


3.37 


00:58:21.87 


+38:00:13.9 


-344.7 


4.98 


-1.04 


-0.69 


22.87 


21.06 


3.06 


00:58:22.02 


+38:04:05.9 


-324.3 


20.49 


-0.84 


2.1 


23.39 


21.84 


1.68 


00:58:23.15 


+37:58:40.9 


-159.6 


8.61 


-2.25 


-1.27 


22.26 


21.1 


4.59 


00:58:24.32 


+38:04:29.9 


-284.3 


2.6 




-3.06 


20.43 


19.71 


2.28 


00:58:25.71 


+38:03:12.2 


-348.1 


11.12 


-0.65 


-0.1 


23.74 


21.53 


2.02 


00:58:29.49 


+38:00:03.6 


-360.7 


19.15 


-0.82 


-0.2 
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-1.09 
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-2.22 
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6.33 
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-361.7 


5.72 


-0.86 


-1.55 


23.38 


21.83 


9.48 



a Distance from the center of EC4 in arcmin. At the distance of EC4, 13 kpc = ldeg. 
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Table 5. Properties of candidate M31 halo stars in stream C2 region (field F36). 
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Table 6. Properties of candidate M31 halo stars in stream B region (field M6). 



a (J2000) 


5 (J2000) 


v r ( km s 1 ) 


Kerr ( km S 1 ) 


[Fc/H] phot 


[Fc/H] spec 


V-mag 


I- mag 


01:08:31.0 


37:30:21.6 


-335.4 


10.1 
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22.51 
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-2.25 


24.38 
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01:09:43.1 


37:41:33.4 


-152.8 


9.844 


-0.58 


-2.1 


23.15 


21.15 
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Table 7. Properties of candidate M31 halo stars in stream A region (field M8). 



a (J2000) 


5 (J2000) 


v r ( km s 1 ) 


Verr ( km S 1 ) 


[Fe/H] phot 


[Fc/H] spec 


V-mag 


I- mag 
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5.416 


-1.41 


-1.55 


22.16 


20.7 



